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Abstract 

I construct the ground state, up to first nonperturbative order, 
of tlie stocfiastic stabilization of the zero dimensional matrix model 
which defines 2D Quantum Gravity. It is given by the linear combina- 
tion of a perturbative wave function and a nonperturbative one. The 
nonperturbative behaviour which arise from the stabilized model and 
from the string equation are similar. I show the modification of the 
loop equation by nonperturbative contribution. 
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1. Introduction 

The zero- dimensional one matrix model leads to a nonperturbative def- 
inition of simplicial 2D quantum gravity coupled to conformal matter with 
c < llH, I, In the continuum limit, the perturbative topological expansion 
in the genus can be evaluated by solving some differential equation: the string 
equation]^]. The Schwinger- Dyson equations of the matrix model defines the 
loop equation of the 2D quantum gravity 0. The matrix model defines also 
the nonperturbative behaviour up to the ambiguity of the boundary con- 
dition of the string equation. Unfortunately, in the case of pure gravity, 
the nonperturbative real solutions of the string equation for every boundary 
condition, can not be solutions of the loop equation^]. 

The stochastic stabilization of the zero-dimensional matrix model is de- 
fined by a one-dimensional matrix model with the same perturbative expan- 
sion in and perhaps it provides a consistent definition of nonperturba- 
tive effects[^, ^. In reference^ the scaling part of the stabilized hamiltoninan 
is solved numerically. In this paper I propose an analytic approach in order 
to show the relationship with the nonperturbative behaviour of the original 
matrix model. 

In section two, I review some results about the nonperturbative effects of 
the matrix model in order to compare it with the stabilized model. 

In section three, the stochastic stabilization of the matrix model is intro- 
duced, the relationship between the eigenvalue density of the matrix model 
and the fermionic density of the stabilized model, and between the Schwinger- 
Dyson equation of the matrix model and the Ward identities of the stabilized 
model are showed. 
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In section four, I perform the calculation of the nonperturbative contri- 
bution in the case of pure gravity. This calculation is performed using the 
old WKB method. Finally, some conclusions are presented. 
2. Nonperturbative 2D quantum gravity 
In this section I review some results of references |P, |10|. 



The nonperturbative effects of the zero- dimensional matrix model are 
given by the tunnelling of one eigenvalue outside the main well of the po- 
tential. The potential of pure gravity is unbounded from below and the 
tunnelling takes place between the well and the unbounded region. 

In the large N limit, the free energy is 

In Z = N^ (^j dXdfip{X)p{fi) In I A - /i I - y dXp{X)W{X)J (1) 

where p(A) is the eigenvalue density and W{X) is the potential. The variation 
of the free energy when one eigenvalue is moved outside the support of p(A) 
is 

r(A/) =6\nZ = N^ J dX6p{X) (2 J dpp{p) \n\X-p \ -W{X)^ (2) 

where Sp{X) = N^^{6{X — A/) — 6{X — Aj)}. One can interpret r(A/) as the 
effective potential for one eigenvalue in the background created by the — 1 
remaining eigenvalues [0]. The formula can be written p, |TT| 

r(A/) = N dX (2 J dpy^ - W'{X)^ 



-2N dAf/(A)VA2 - a? (3) 



where a is the cut of the support of p(A) and U{X) is a polynomial given by 
the saddle point aproximation. The maxima of the effective potential r(Aj) 



are given by the zeros of U{X). One may consider the instanton configuration 
where one eigenvalue sits at the top of r(A/). Taken into account that the 
eigenvalue density is ITT 

p(A) = U{X)Va^ - A2 , A G (-a, a) (4) 

one can write the instantonic contribution as an integral of the imaginary 
part of the analytic continuation of the eigenvalue density 

Tmatr^. = -2N [' dX^piX)) (5) 
J a 

where a is the cut of the support of p{X) and b is the top of the effective 
potential r(A/). In the double scaling limit this contribution is finite and is 
the source of all the troubles of pure quantum gravity. This nonperturbative 
contribution arise also from the string equation 0. 

The instantonic action (Q) can be interperted as an effective potential as 
follows: a source coupled to an eigenvalue is added to the original potential 

N N 

5^Vr(A,) -^5:iy(A,)- JA;v (6) 

i=l i=l 

At leading order the source term only affect the eigenvalue Aat, so that the 
eigenvalue density is now 

p(A) = p(A) + 1{<5(A - A^) - 6{X - A^)} (7) 

where Aat is the eigenvalue of the saddle point configuration in the presence 
of the source and X^ is the eigenvalue without the source. In the large N limit 
the classical field (Atv) is given by the eigenvalue X^. Hence, the effective 
potential 

r((A;v)) = lnZ(J)-(JA^) (8) 
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becomes the instantonic action when the classical field and A/ are iden- 



tified. 



It is well know that SU{N) symmetry is broken if there exist non trivial 
solution of 



Henceforth, the instantonic configurations with an eigenvalue at the top of 
the effective potential break the SU{N) symmetry. 

Therefore, SU{N) symmetry is restored like follow: there are vacua \i) 
with an eigenvalue Aj at the top of the effective potential. SU(N) vacua are 
given by linear combination of \i) which must be invariant under permutation 
of the eigenvalues. There are only one linear combination and does not 
exist arbitrary parameter like the 6 parameter of gauge theories. 

3. Stochastic stabilization of the zero-dimensional matrix model 

Stochastic stabilization ^ provides a nonperturbative definition of 2D 
quantum gravity, while reproducing the perturbative expansion in powers of 
of the original matrix model and, therefore, the topological expansion 
of the discretized quantum gravity. 

The stochastic stabilization introduces a positive definite hamiltonian 



this hamiltonian is well defined for all values of and coupling constants. 
The zero mode of the stabilized theory is 



its norm being the partition function of the original matrix model. Hence, 
the matrix model is well defined if and only if this zero energy state is a 
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(10) 




(11) 



normalizable state. When this is the case, corresponding obsevables in both 
theories coincide 

{Q)stab = j Q\-^\^ d<^ = ^ j gexp{-W($)}rf<l> = (g)„„t„. (12) 

In the A„ variables, where {A„} are the eigenvalues of the matrix $, the 
zero energy state is the Slater determinant 

N { N ^ 1 

^({An}) = n(A. - A,) exp -- W{Xn) (13) 

i<k \ ^ n=l J 

and, hence, the stabilized model is a Fermi gas. 



In the planar limit, the condition ( [T^ ) becomes 



/ dW"' Pmatrix = / dXX^ p stabilized (14) 

where Pmatrix and pstabiUzed are the density of eigenvalues of the zero dimen- 
sional matrix model and the fermionic density of the stabilized model. And 
taking into account that both densities are the square root of a polynomial 
which is zero at the points a and h it is not difficult to prove that both 
densities must be equal0. 

The fermionic density in the large limit for the stabilized hamiltonian 

is 
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p{X) = -p{EF-V) (15) 

where V is the stabilized potential^ and Ep is the Fermi level. From the 
equality between ([TsD and (^, Bp — V must has a simple zero for A = a and 
double zeros corresponding to zeros of f/(A), in particular, a double zero for 
X = b. Therefore, in the simplest case of a cuartic potential |jl2|, the stabilized 



^in the large N limit, V is a one particle potential given by the Hartree aproximation 



5 



model is given by a Fermi gas placed at the main well of the stabilized 
potential, where the classical turning point of the Fermi level is given by the 
cut of the eigenvalue density of the matrix model: a. The stabilized potential 
has also a local minimun 6, which is reached by the Fermi level, and it is given 
by the top of the effective potential of the matrix model. Hence, the main well 
of the stabilized potential defines the perturbative expansion of the matrix 
model and I will show how the local minimun defines the nonperturbative 
effects. 

The stabilized hamiltonian ( p!0D is the supersymmetric hamiltonian of 
reference restricted to the bosonic sector. The Schwinger-Dyson equation 
of the zero-dimensional matrix model becomes the Ward identities of the 
supersymmetric matrix model. The Ward identities of the stabilized hamil- 
tonian are given by the condition ||13|| 



where gn are the coupling constants of the zero- dimensional matrix model 
and E is the ground energy of the stabilized hamiltonian. In appendix 2, I 
extract from this condition the Virasoro constraints. 

4. Nonperturbative stochastic pure gravity 

Let us consider the matrix model 



and $ is an dimensional hermitian matrix. This model corresponds to a 




(17) 



where 



(18) 



discretized formulation of pure quantum gravity 0. 
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The Fokker-Planck hamiltonian of the stabihzed model is the sum of 
one particle hamiltonians ^ 

1 



and an interacting term 



(19) 



(20) 



where {A„} are the eigenvalues of the matrix $. The interacting term is 
subleading in the 1/N expansion and I discard them. The one particle stabi- 
lized potential (|19D has a main well and a local minimun below some critical 
coupling constant and only one well above it. 

The expansion is equivalent to the WKB aproximation, and the 
quantization condition is0 

r I 1 

- / d\J2{E^ -v) + --- = n + - (21) 

The Fermi level E^^i in the large limit is the value of the one particle 
potential (0) in its local minimun 1 12]. The first correction to the Fermi 
level is negative (appendix 3). Henceforth, in the perturbative expansion 



1/A^, the quantization condition (21) restricted to the main well gives the 
correct perturbative quantization of energy levels. 

However, the perturbative WKB wave function has a singularity in the 
local minimun. It is well know that the WKB aproximation break down near 
the turning points where the variation of the local wavelength is not small 

" ■ ■ (22) 




'2{E-V)^ 

At those points one needs the solution of the Schrodinger equation in order 
to give the connection formulas and the quantization condition. 
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The local wavelength near the local minimun of the stabilized potential 
(HP is 

A(a:) ~ I (23) 



where eih is the perturbative shift of the Fermi level and ei > because the 
Fermi level is placed below the local minimun. The origin of coordinates is 
placed at the local minimun and x = X — b. 



The variation of the local wavelength (pB]) is order one when x is order 
y/h. Henceforth, I perform the change of variables 

and the Schrodinger equation at first order in h becomes 

^^-(ei + CV)^ = (25) 

In appendix 1, I show the solutions of this equation and the connection 
formulas in the local minimun which arise from them. 

The more general wave function in the main well of the stabilized potential 
is the oscillatory solution 

*<" = i5^^r?p~K'^r''-v/^(^ + *). (26) 

Because the potential is even, one can construct a wave function with definite 
parity. For the even wave function 5 = and 5 = | for the odd wave function. 

The energy level may be written as the sum of a perturbative and a 
nonperturbative term : E = Cn + Cn- Hence, the wave function becomes 
^(A) = ^P(A) + *"P(A), where 
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where 

r? = ^ - iV £dx^2{er,-V)-S, (28) 

a is the classical turning point and a;„ is the classical frecuency. The wave 
amplitude of ^"^(A) is suppressed by the nonperturbative term e„, hence 
^^(A) is the perturbative part of the wave function. 

Beyond the local minimun the wave function must be the fall off expo- 
nential in the infinity 

where b is the local minimun. 

The connection formula at the local minimun gives the wave function 
between the local minimun and the main well (appendix 1) 



where Xq is an arbitrary constant. This wave function may be written as 
^(A) = ^+(A) + ^_(A) where 

^+(A) A 



Bfi{ei,Xo) {2{V-E))4 



r exp I + £ dx^2{V-E)^ 
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b 



A = exp^-NJ dx^2{V-E)j (31) 

where \1/+(A) is suppresed by a nonperturbative term, hence the perturbative 
part of the wave function is \1/_(A). 

The usual matching condition in the classical turning point a between 
the perturbative wave functions \E'^(A) and ^I'-(A) gives the perturbative 
quatization condition (^Tj) and the following relation between the arbitrary 
constants A and B 

A = 2Bf2{ei,Xo)expi^N dx^2{V - E)"^ (32) 

And the matching condition between \E'"^(A) and ^E'+(A) gives the follow- 
ing value for the nonperturbative part of the energy level 
1 /(ei,xo) 



J a 



b dx 



■ exp |-2A^ dx^2{y - en) | (33) 



where, /(ei,Xo) ~ Xq for xq small (appendix 1) . This nonperturbative 
corrections arise only for levels near the local minimun. 
Therefore, the total energy of the ground state is 

^ = E en (34) 

n=0 

where e„ is nonzero only for levels near the Fermi energy where 

en = EF + 0{l/N) (35) 
hence, in the large limit, ( p4|) is approached by 



E = ^ exp {-2N dXJ2{V - Ep)] V , ^ (36) 
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The nonperturbative contribution to the observables can be written as 

K{g)exp{T stabilized), (37) 
where the instantonic action is given by 

TstaMHzed = -2N f d\%p{\)) (38) 

where p(A) is the analytic continuation of the fermionic density of the stabi- 
hzed model which is equal to the eigenvalue density of the zero-dimensional 
matrix model. Hence, the instantonic action of the matrix model (|^) and of 
the stabilized model (|38|) are equal. (|38|) also agree with the instantonic ac- 
tion of reference^], which is given by a sucession of instanton-antiinstanton 
starting and ending at the main well. 

The constant Xq is analogous to the arbitrary constant of the nonpertur- 
bative contribution which arise from the string equation 0. However, in the 
stabilized model there is not a relationship between xq and boundary condi- 
tions of some string equation, in fact the ambiguity arise form the connection 
formulas (appendix 1). In reference]^ the stabilized model is solved numer- 
ically and the solution is unique, so the ambiguity which I have found must 
be a defect of the approach, but the arbitrary constant Xq must be different 
from zero because otherwise the WKB wave function is not a solution of the 
Schrodinger equation near the local minimun. 

The important point is that the ground state energy is greater than zero 
by a nonperturbative correction. From reference []13| and (|3^) the first loop 
equation becomes in the double scaling limit 

V {W{L)) - dJ{W{J)W{L - J)) ~ f{z) exp |-^v^^l | (1 + . • .) 

(39) 
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hence, the nonperturbative behaviour of this new loop equation is compatible 
with the string equation of the zero dimensional matrix model. But, in the 
stabilized model the string equation and the KdV flow are lost because the 
Virasoro constraints also changes by nonperturbative corrections. Is an open 
problem what replaces the old string equation and the KdV flow. 
5. Conclusions 

I have found the first nonperturbative correction to the ground state of 
the stochastic stabilization of zero dimensional matrix model. The nonper- 
turbative behaviour of the original matrix model and its stabilization are 
similar. 

The loop equation of the stabilized model and the zero-dimensional ma- 
trix model differ by a nonperturbative term (|39|), which is similar to the 
nonperturbative behaviour of the string equation. Therefore, I expect that 
the stabilized model gives a consistent definition of the simplicial 2D quan- 
tum gravity. 

From the point of view of the original matrix model, the nonperturbative 
effects arise from the tunnelling between the well of the potential to the un- 
bounded region. Hence, the perturbative vacuum state is a metastable state 
and it decay to an ill defined vacuum. But, from the point of view of the 
stabilized model there is only one perturbative vacuum state and there is 
not an ill defined vacuum [|l^, ^. In ordinary quantum mechanics the non- 
perturbative effects arise from tunnelling between two or more perturbative 
vacua, the true nonperturbative vacuum is given by some linear combination 
of perturbative vacua, and the physical interpretation is that one pertur- 
bative vacuum decay by quantum tunnelling into the others. However, the 
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ground state of the stablized model is given by the sum of a perturbative 
part, which is the perturbative wave function around the main well, and a 
nonperturbative term 



^(A) = ^P(A) + ^"P(A). 



(40) 



The perturbative wave function is given by an oscillatory solution at the 
main well and a fall off exponential at the classical forbidden region 

A 



^P(A) 
^P(A) 



- cos [n r dxJ2(E - V) -- + ri\x<a 

exp I -A^ dxJ2(V-E) \ \ >a (41) 

Therefore, the particle density at the main well is given by 

p{\, E) = ^ cos' (n r dyp{E-V{y))) . (42) 

In the large N limit the cos' must be replace by 1/2 and the particle density 



is given by the WKB formula p4 



p{\E) 



p{E-V{\)) 



(43) 



But, in the double scaling limit (|42|) becomes 



p(x,e) 



1 



: COS" I ^ 



dy\/2{e -v{y)) 



(44) 



^2(6 -vix)) 

where e and v are the scaling part of the energy and the potential, and h is 



the scaling coupling constant |n 



N'ig-gc 



(45) 
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where Qc is the critical couphng constant. The argument of the cos^ is a 
function of the scahng constant h, and now the particle density becomes an 
oscillatory function if h is finite. The amplitude is an increasing function 
of the position and go to infinity when the position approach the classical 
turning point because the WKB approach break down at the turning points. 
In fact the exact particle density is computed numerically in|^ and it is given 
by an oscillatory solution with an increasing amplitude when the position 
approach the turning point, but, of course, is finite at the turning point. 
Therefore, the oscillations of the particle density, which are nonperturbative 
in nature, appears at the first WKB aproximation of the wave function. 

The nonperturnative part is given by the same oscillatory solution at the 
main well with a phase sifted by 7r/2, and an increasing exponential between 
the cut a and the local minimun b, and beyond the local minimun there is 
only the perturbative wave function. 

^"^'(A) = sin In r dxJ2(E -V)-- + r]\x<a 

^"^A) = -^^^^^^^exp[+Njyx,/W^)}a<X<b (46) 

The amplitude of the nonperturbative wave function is suppressed by a non- 
perturbative term 

A ~ exp |-2A^ dx^2{y - E)^A (47) 

In reference the particle density beyond the local minimun h decrease 
very fast. In the double scaling limit the perturbative wave function (p]) 
decrease as 

^'~exp|-^y" ^2{y - e)| ~ exp (48) 
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Between the cut a and the local minimun b the increasing exponential of 
the nonperturbative wave function gives some contribution to the particle 
density. Hence, the particle density decrease more slowly between a and b. 
However, the nonperturbative wave function is supressed by a nonpertur- 
bative term, and the particle density also decrease exponentially between a 
and b. This is in agreement with]^ where the particle density decrease more 
slowly between the cut a and the local minimun b than beyond the local 
minimun. 

The above discussion suggest that there are only one perturbative vac- 
uum, given by and almost one nonperturbative vacuum, given by 
So the perturbative vacuum decay into a nonperturbative vacuum. The non- 
perturbative vacuum break down the symmetries of the perturbative one, so 
the Ward identities of the model (loop equation and Virasoro constraints) 
must change. Is an open question what kind of symmetries replaces the old 
one. 
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Appendix 1 

In a neighbourhood of size around the local minimun the Schrodinger 
equation is given by 

0-(ei + CV)vl> = O (49) 

where ei is the first perturbative correction to the Fermi level and is positive 
because the perturbative Fermi level is placed below the local minimun. I 
have placed the origin of coordinates at the local minimun of the potential 
and X = X — b. 

Solutions of (|i9| ) are given by linear combinations of the even and odd 
functions 

^o..(x) = v^a:e-^^^>(^ + ^ I ^ I Cx^) (50) 

where F{a \ c \ z) is the Hipergeometric degenerate functions and is the 
solution of|T^ 

zF" + (c - z)F' -aF = (51) 

The asymptotic behaviour of (0) can not match the WKB wave function. 
However for large values of x, ei is small in comparison with the potential, 
and one can approach (|50| ) by Bessel functions of order ^ at zero order in 
ei. The asymptotic behaviour of Bessel functions match the WKB wave 
functions. 

Hence, I construct a wave function, with the correct asymptotic behaviour 
and given by (0) around the local minimun, as follows. 

The real axis is divided into three regions: the neighbourhood of the local 
minimun, where the wave function is given by Hipergeometric degenerate 
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functions 

^1 = a(^e^en + ^odd) + Pi^even " ^odd)- (52) 

An intermediate region, where the wave function is given by Bessel functions 



'—X 



vl>+(x) 



X 



= A2 

+ B2 

= 

+ 5i 



, Cx , , Cx? , 

4V 2 ^ — V 2 • 



4^ 2 



— V 2 ' 



(53) 



The continuity of the wave function and its derivative is required in the 
boundary xq of the above regions 

^i(xo) = ^^(xo) 



^l(-Xo) 

dx 



x=xo 



'd^ 

dx 



\ dx 

(d^ 

\ dx 



x=xo 



(54) 

x=—xo \ / X=—Xq 

And the asymptotic region, where the Bessel functions becomes WKB wave 
functions. The asymptotic behaviour of Bessel functions are 

2 , 



Ii{z)+I_i{z) 

4 ^ ' 4 ^ ' 



n(z) 



V2^ 

-2 . TT _ 

—j= sm(-)e 



(55) 



The set of equations ( [5^ ) gives the connections formulas between the 
asymptotic WKB wave functions. 
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The wave function must be the fall off exponential for x ^ 0. Hence, Ai = 
0, and Bi = B is an arbitrary constant, which can be fixed by normalization. 
For a; ^ the wave function is 



where /i and /2 are given by the condition (@). 

I have constructed a wave function which is given by the usual WKB 
wave function in the asymptotic region, and in the interval (— a;o,a;o), by the 
exact solution of the aproximate Schrodinger equation (^). However, this 
wave function is not unique because the point xq is arbitrary. 

In the limit of Xq small, the coefficients of the wave function ( ^6]) becomes 

1 

/i(^o) = ^0 (r(5/4))2 \/7 + ^^^o) 

/2M = l + ^o^m^/f + (57) 

If one set Xq = 0, the wave function for x <^ is given only by the increas- 
ing exponential. This connection condition also arise from the Schrodinger 
equation (^) if one set ei = 0. Hence, the decreasing exponential can be 
interpreted as the first correction to the approach ei = 0. 
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Appendix 2 

In this appendix I show the proof of the Virasoro constraints in the 
stochastic stabihzation of 2D quantum gravity. The proof of the Loop equa- 
tion is given in reference [[T^j . 

Let be the matrix potential W 

W = Y^ '^Tr^^ (58) 

n=0 

the potential of the stabilized theory becomes now 

-| I oo oo 1 CO n— 2 I 

^ = ^ E E ^?n^7™Tr$"+'"-2 - T7 E E Tr^'Tr^'-'-' (59) 

^ [n=Om=0 n=0 p=0 J 

and from the Hellmann-Feynmam theorem 

arp -| I oo 1 k—2 I 

— = - 2 ^ gr.{Tr<^-^'-') - -J2{Tr<!>^Tr<!^'^-^-') (60) 

(ygk ^ I n=o p=o I 



If the ground energy is zero, (|60D becomes the usual discrete Virasoro 
constraints iTBI 



o d 1 d"^ 

J2 9n{k + n- 2)- ^ Y^p^k-p- 2)-— I Z = 

,^0 dgn+k-2 4iV2^o dgpdgk-p^2' 



(61) 
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Appendix 3 

In this appendix I perform the calculation of the first correction to the 
Fermi level in the stabilized potential of the cuartic matrix model and I will 
show that the Fermi level is placed below the local minimun. 

The perturbative quantization condition for the Fermi level is given by 



^ r dxp{E,-V)-l- r dx +0{llN) = N-\ (62) 



where the extra term arise from the Hartree aproximation to the interacting 
term (^Ol). The Fock term is subleading in the quantization condition (|6^). 
The first correction to the Fermi level is given by 

where uj is the classical frecuency which is positive. 
Following reference [p^] E^P becomes 

Ef ^ I a" ( / u \\ 



-1 + == {71-2 arctan , ) (64) 

and is easy to see that for 6^ > > 0, e'P is negative and the Fermi level 
is placed below the local minimun. 
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